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1. Introduction 



In recent years, much attention has been paid to the short-time critical dy- 
namics. The short-time phenomena arise at times just after a microscopic 
time scale t m i C needed by the system to remember only the macroscopic con- 
dition and to forget all specific microscopic details. The corresponding time 
regime is also called critical initial slip in order to distinguish it from the 
uninteresting microscopic time interval between zero and t mic . Since the pi- 
oneering analytical study of universal short-time scaling has been found 
in various models ( See || and ||). When the system is quenched from a 
high temperature Tj to the critical temperature T c <C Tj the order parameter 
shows in the short-time regime a power law increase m(t) ~ t 9 with a new 
universal critical exponent 9' . 

The short-time dynamics has been thoroughly investigated for models 
with short-range interaction (SRI). Since the critical equilibrium properties 
are modified by the presence of long-range interactions (LRI) it may be in- 
teresting to know how the short-time critical behaviour depends upon the 
interaction range. 

The statistical mechanics of LRI has a long history. Already in the 60's 
Jovce 0] studied thermodynamic properties of the static spherical model with 
long-range ferromagnetic interaction between the spins. The static critical 
exponents for LRI have been computed for the n- vector model by use of the 
renormalisation group approach [4-8] and the 1/n-expansion techniques [[K] 



There are also Monte Carlo simulations for the one- dimensional static model 
The dynamic properties of the LRI in the long-time regime have also been 



studied as early as the 70's. Suzuki et al. |L2[ extended the dynamic theory 



developed by Halperin et al. [|13j to investigate an exponent which describes 
the critical slowing down in the n-vector model with LRI for T > T c , with 
equilibrium initial conditions. Folk and Moser studied a three-dimensional 
dynamical model for liquids and demonstrated that the critical dynamics was 



affected by LRI [14]. The kinetic spherical model showed that the short-time 



critical exponents were modified by LRI [|i5] . 

In this paper, we study the short-time critical behaviour of the dynamic 
Ginzburg-Landau model with long-range exchange interaction. In equilib- 
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rium at temperature T the 0(n) symmetric Hamiltonian is given by 



|(v s ) 2 + ^(v* s ) 2 + ^ 2 + f( s 2 ) 2 } (1) 



where s = (s a ) are n-component spin fields, r is proportional to the reduced 
temperature 1 — T/T c and g is the coupling constant. The SRI model cor- 
responds to a = 1 and 6 = 0, whereas for the pure LRI model a < 2, a = 
and 6=1. Since the case 0<cr<<i/2is covered by a mean- field theoretic 
description, and since for o > 2 and d > 2 the model (1) belongs to the same 
universality class as the SRI model, we will restrict ourselves in the present 
paper to the range d/2 < o < min(2, d). 

The dynamics to be discussed here is controlled by the Langevin equation 

where A is the kinetic coefficient. The random forces £ = (£ a ) are assumed 
to be Gaussian distributed 

< £ a ( x , t) >= ; < C{x, t)f{x', t') >= 2X5 a/3 5{x - x')5(t - t') . 

As mentioned above, the initial state is prepared (macroscopically) at 
some very high temperature Tj. One assumes that the initial condition 
s (x) = s(x, 0) has also a Gaussian distribution P[s Q ] oc exp(— iP[s ]) where 



H l [s ] = J d d x^[s (x) - m (x)f 



To is proportional to 1 — Tj/T and mo(i) is the (spatially varying) initial order 
parameter. Being away from criticality (Tj ^> T c ), the initial correlation 
function will be short-ranged. Since tq ~ \f (where \i is a renormalisation 
momentum scale), the physically interesting fixed point is Tq = +oo, which 
corresponds to a sharply prepared initial state with initial order mo and zero 
correlation length. 

Introducing a (purely imaginary) response field s(x, t) |T7|, [T8j , the gener- 



ating functional for all connected correlation and response functions is given 
by 

W[h, h} = \nj V(is, s) exp j -C[s, s] - H%s ] + J dt J d d x(hs + hs) 1 (2) 
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where 



OO s 

C[s, s] = J dt J d d x\s h + A (t - a y 2 +&(-y 2 )^) s + 



^9 2 
~6 SS 



Xs 2 



step function B(t = 0) = 0. Then the contribution 
jC[s, s] arising from the functional determinant det 



vanishes [IIT 



(3) 

Here we have used a pre-point discretisation with respect to time so that the 

proportional to 6(0)) to 

5s(x, t) 

It is believed that the singularity of the temporal correlation is essential 
to the short-time scaling and the scaling can emerge in the early stage of the 
evolution even though all spatial correlations are still short-ranged. 

The system is now rapidly quenched to a temperature T ~ T c . The order 
parameter will undergo a relaxation process displaying an initial increase. 
As long as the correlations are short-ranged and the spatial dimension d is 
smaller than the critical dimension d c , the order parameter follows a mean- 
field ordering process because the mean-field critical temperature Tj mf ^ is 
larger than the actual critical temperature T c . This ordering causes an am- 
plification of the initial order parameter. For d > d c mean-field theory applies 
and there is no critical increase. 

For the SRI models d c = 4. The longer is the interaction range, the 
stronger the suppression of the fluctuations and hence the critical dimension 
of the LRI model is smaller. Indeed, it turns out that d c = 2a. Also one would 
expect that the critical initial increase should be weaker as the interaction 
range becomes longer. 

Since the short-range exchange interaction is irrelevant for d/2 < a < 
a s = 2 — f] sr where 7] sr is the Fisher exponent at the SRI fixed point, one can 
consider only pure LRI. We apply the e-expansion theory to the LRI model 
in this regime with e = 2a — d. The critical initial order increase appears in 
the LRI model for 1 < d < d c . The scaling behaviour of the critical initial 
slip is governed by the exponents 9 and 9' . They are computed as functions 
of d and a. 

For a close to (but larger than) d/2 the quantity e is very small and the 
numerical values of the exponents are accurate when computed to order e 2 . 

However, when the interaction range is not very long, the situation be- 
comes more complicated, due to a subtle competition between the SRI and 
the LRI fixed points [6-9]. Honkonen || computed the /3-function of the 
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renormalisation group for the pure LRI model at three-loops and found that 
the infrared LRI fixed point becomes unstable for o = a s . In the pure LRI 
model the exchange interaction term is not renormalised, so that the anoma- 
lous dimension of the field s(x, t) is zero, whereas in the SRI model the field 
carries some anomalous dimension 7. Taking the limit o — > 2 the expres- 
sions for the anomalous dimension (and for other critical exponents) do not 
coincide. However, as first shown by Sak || to the leading non-trivial or- 
der and later by Honkonen and Nalimov |7]] to all order in e' = 4 — d, the 
anomalous dimension 7 and the other exponents are continuous functions of 
the parameter a. This means that the scaling regime of the LRI model is 
valid only for a < a s , whereas for a > a s the scaling behaviour is described 
by the SRI model. At the borderline value a = a s the two descriptions yield 
equal values for the critical exponents. Let us conclude here by remarking 
that these last results were obtained solely for static models. 

The paper is organised as follows: In section 2, the LRI model with a < a s 
is studied by the e expansion method. The scaling behaviour of the order 
parameter, correlation and response functions, as well as the corresponding 
critical initial slip exponents, are obtained. Section 3 contains conclusions 
and discussions. 



2. The short-time scalings and exponents 

Since the SRI is irrelevant for a < cr s , in this section we take a = and 6=1 
in®- 

For g = 0, the generating functional (2) becomes Gaussian and can be 
easily evaluated in momentum space. One must take into account the initial 
condition, by imposing the following boundary conditions: 

s(x, 00) = sq(x) = mo(:r ) + r _1 s(x, 0) . 

The free response function G p (t,t') and the free correlation function 
C p (t,t') are respectively 

G p (t,t') = e(t-t')exp[-\(p a + T)(t-t')}] 
C p (t,t') = C^(t-t')+C^(t,t'), 
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with equilibrium part C^\t — t') and initial part C^'(t,t') defined by 

C£\t-H) 
C®(t,t') 

One sets now a perturbation expansion ordered by the number of loops in 
the Feynman diagrams. It is convenient to consider the Dirichlet boundary 
conditions To = +00 and mo(x) = 0. The general case is recovered by treating 
the parameters r _1 and m (x) as additional perturbations. 

The model (2) with Dirichlet boundary conditions must be renormalised. 
For this purpose notice that the free correlation function simplifies to 

Cf \t, = — ^ {exp[-A(p CT + r)\t- t'\] - exp[-A(p" + r){t + t')}} . 

By integrating over the internal momentum and time coordinates one en- 
counters ultraviolet divergences which can be absorbed through the repa- 
rameterization of a finite number of coupling constants and fields. 

Through dimensional analysis, one can show that the critical dimension 
d c = 2a and hence it is convenient to make an expansion in e = la — d. We 
will adopt the dimensional regularisation with minimal subtraction scheme 
|p0|1 and introduce renormalised quantities through multiplicative factors 

s b = Z\l 2 s, s b = Z± /2 s, X b = (Z s /Z~) l/2 \ , 
n = Z~ x Z t t , g b = K^ l i/Z~ 2 Z u u , 

r ob = (Z 7 /Z s )^ 2 t , s 0b = (Z 7 Z )^ 2 s (4) 

where the subscript b denotes the bare quantity and = 2 1 ~ d n~% [r(<i/2)] _1 . 
Some comments are in order: 

(i) The graphs containing only the equilibrium part of the correlation 
function are associated to 1PI diagrams and can be made finite by the same 
renormalisation factors as the translation-ally invariant theory. 

(ii) In addition there are divergences arising for t + 1' — from the initial 
part of the correlation function. Remarkably enough, such divergences can be 
multiplicatively removed if one associates them with n-point connected Green 



r + p 



-exp[-A(^ + r)|t-t'|] 



T +p a 



exp[-X(p a + T)(t + t')} . 
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functions. A simple dimensional analysis reveals that new re-normalisations 
are required only in two-point functions. Due to the Ward identities: 



s (x) = 



s (x) = 2\s Q (x) 



which hold when inserted in the connected Green functions, one is left with 
a single additional renormalisation constant Z . 

A two-loop calculation gives the following renormalisation constants: 



i; 

1 n+2 R 2 

1 t, — B a u ; 

6e 

n + 8 \(n + 8) 
1 H t, — u + 



Z T = 1 + 



6e 
n + 2 



hn + 22 



6e 



-u + 



36e 2 36e 
(n + 2)(n + 5) n + 2 
24T 



Z = 1 + 



n 



6e 



2 n + 2 

— u H — 

12e 2 



36e 2 
n + 5 



-Da 



It 



u 



(5) 
(6) 

(7) 
(8) 
(9) 



Here we have introduced 



d 2a x 



[1 + x a + (e + x) 



^- 2 x~ a 



(2tt) 2ct 

with e a unit vector in the 2<r-dimensional space, and 

£ ff = ^(l)-2#7/2)+#7) 

with the logarithmic derivative of the gamma function. For the partic- 
ular case a = 2, one has B 2 = |ln(4/3), and D 2 = 1. The calculation of the 
renormalisation constant Zq is reported in Appendix. 

According to the general solution of the renormalisation group equation, 
the renormalised connected Green function of N s-fields, N s-fields, and M 
s -fields at the fixed point u* has the following scaling law: 

G M ~({X, t}, T, Tq 1 , A, U*, /i) = l( d -v+Vs)N/2+(d+a+ri s )N/2+(d+a+ Vs + Vo )M/2 



X G^~({IX, l a+ « U *k}, Tl~ a+ < U *\ TvH a+ « u *\\, U\ U) 



(10) 
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where rj s = j(u*), rjs = j(u*), and 770 = 70 («*)• The Wilson functions 
entering the renormalisation group equations are defined by 



7 = fidp lnZ s \ ; (3 = fid^u] 



> 



7 = /i9„lnZ 7 | ; 



k = //^lnr| ; C = /i<9 M lnA| = 5(7 - 7) ; 70^/^lnZ, 







and are computed perturbatively from Eqs.(5-9). The symbol |o means that 
//-derivatives are calculated at fixed bare parameters. For instance, at the 
two-loop level, the Wilson function 70 (related to the initial order parameter) 
is given by 

n + 2 



7o 



6 



1 + 



-In2-^D„) „ 

a 2 



u . 



(11) 



By solving algebraically the equation j3{u) = one finds the infrared LRI 
fixed point 



6e 



u 



n + 8 



2(5n + 22) 
+ (n + 8) 2 CT 



+ 0(e 



(12) 



and subsequently the values of the Wilson functions at this point. 

In order to identify the critical exponents one can compare the standard 
scaling form of the two-point correlation function 



r\-(d-2+r)) 



f 



/"V"» ry*' rv* ry>' rv* iy' 

\Ju tXj JU JLj JU J> 



t ,l/z 



to the Eq.(10) in which we have set N = 2, N = M = and Ix = 1. Here 



In this way we find the LRI critical exponents to second order in e |5], [12 
2 - a + 77. 



V 



a 



z = a + C« 

l/l/ EE CT — k(m*) = (J 



2-c- ; 
6(n + 2 



(n + 8 
n + 2 
n + 8 



2 ^ e 5 



, 7n + 20 ^ 

1 + , . ^o AxC 



n + 8 



€ . 



Notice that the anomalous dimensions of s and s are i] s = 2 — o + 77 and 
?7j = 77 + 2(z — a) respectively, whereas that of the initial order parameter r) 
is given by Eq.(ll) at the fixed point (12). 
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Employing a short-time expansion of the fields s(x, t) and s(x, t), as done 
in |l|], one can derive the following behaviour of the full response and corre- 
lation functions for t > but t' — > 0: 



G( P , t, H) 
C(p, t, f) 



P 



P 



-2+ri+z 



-2+77 



fc(p£,p z t) • 



(13) 



Here we defined the initial slip exponent 9 and computed it to second order 
in e 



e(n + 2) 



2z 2a n + 8 



1 + 



7n + 20 

(n + 8) 2 



12 In 2 



a(n + 8) 



Let us discuss now the scaling form of the order parameter which relaxes 
from a non-zero initial value mo to zero. As mentioned above we can consider 
mo(i) an additional time independent source coupled to the initial response 
field so(x). Owing to the renormalisation of the initial order parameter 



m ob (x 



(ZoZ~ 



-1/2 



m {x) 



no new renormalisation is required for the time dependent order parameter 
m(x,t) =< s(x, t) > \j i=h=Q . By taking a homogeneous source m (x) = m , 

-foo, we obtain the power law 



but keeping still Tq 

m(t) = m t e 'f m (m t 
where the exponent 6' is defined by 



!/ 1 d— 2 + rj 



2z ; ft vz 



(14) 



Vs + Vs + Vo 



2z 



To second order in e it has the value 



e 1 



em 



2a(n + 8) 



1 + 



7n + 20 

(n + 8) 2 



12(^2-0"^) 
oTn + 8) 



As first indicated in [l[] for the SRI model, the critical exponents 9 and 9' are 
related by 9 = 9 + (2 — z — rj)/z. 

The function f m {x, y) appearing in (14) has a universal behaviour at the 
critical point r = 0: / m (0, 0) is finite; while for x —>■ 00, f m {x, 0) behaves like 
~ 1/x. 
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3. Discussions and conclusions 



When the long range interactions are dominant, e is small enough and the 
calculated values of 9' and 9 for physical dimensions are numerically reliable. 
For instance, we list in Table [I] the values corresponding to e = 0.1 for n = 1 
and d = 1, 2, 3. 





d=l, a = 0.55 


d = 2, a= 1.05 


d = 3, a = 1.55 


9' 


0.0383 


0.0180 


0.0117 


9 


0.0408 


0.0187 


0.0120 



Table 1: The values of 9', 9 to e = 0.1 for n = 1 and d = 1, 2, 3. 



In the following, we will focus the discussion on 9'. Let us first notice that 
both the response and the correlation functions measure the fluctuations of 
the spin- fields. Since 9 and 9' are positive, one expects, according to Eqs.(13), 
an initial increase of the fluctuations. Of course, the increase depends upon 
cr and d. 

In Figure 1 the exponent 9' is plotted versus d for a — 1/2, 1, 3/2 and 2 
respectively and n — 1. The value a = 2 corresponds to the SRI model. At 
fixed a, the exponent 9' decreases when d increases, because fluctuations are 
reduced as the dimension gets larger. At the critical dimension d c = 2a the 
value of 9' becomes equal to zero. Here other scaling laws would replace the 
power law. 

Figure 2 shows that 9' for d = 2 and small a monotonously increases with 
n. For larger a it first reaches a peak and then it decreases toward n — > oo. 
For other values of the spatial dimension d the pictures are similar. The 
increase of & can be easily understood as more internal degrees of freedom 
(larger n) help the fluctuations increase. Hence the critical behaviour is 
smooth in n and can be studied in an 1/n-expansion. But for large a one 



can reach the opposite effect |2T|, the fluctuations decrease when n exceeds 
some threshold value. 

In Figure 3 the exponent 9 is plotted versus a for n = 1 and d = 1, 2, 3. At 
fixed d the exponent 9' decreases when a decreases, because the fluctuations 
are more suppressed by interactions of longer range (a smaller). In one 
dimension, there is no SRI fixed point hence only the curve controlled by the 
LRI fixed point is observed. 
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All the previous considerations were of qualitative nature since they do 
not take into account the interaction of coupling constant g. Of special 
interest in this respect is the LRI fixed point (12). At a = a s = 2 — i] sr 

(where n sr = — — — -r-e' 2 and e' = 4 — d) and fixed d we have 



2(n + 8) 2 



n + 2 



4(n + 8) 



e' 2 - 



6e' 



u = 



6' 



n + 8 
e'(n + 2 
4(n + 8 



1 + 



3(3n + 14)e' 



1 + + 



(n + 8) 2 
6e' 



it 



Si? i 



n + 3 , 3 
, h hi - 

n + 8 Vn + 8 2 



= 9' 



SR ■ 



(15) 



Here the subscript SR means short-range regime. 

In order to explore the limitation of o — > 2 of the LRI, we make a double 
expansion in e and a = 1 — a/2 with a of the order e or smaller. The infrared 
fixed point to order e 2 is located at 



6e 



u 



wlr 



n + 8 



1 + 



n + 8) J 



3(3n + 14) + (n + 2) 



a 



a + e 



(16) 



Here the subscript w/r (weak-long-range) means that a is at most of order 
e. The critical initial slip exponent in the weak-long-range limit can be also 
computed to this order: 



& 



wlr 



e(n + 2) 
4(n + 8) 



1 + a + 



n + 8 



6(^ + ln^ 

n + 8 2 



+ 



n + 2 a 
n + 8 a + e 



(17) 



Since a is actually of order e 2 it can be set to zero in (16) and (17), and 
taking into account that e = e' + 0(e 2 ), one gets 



u lwr 
v Iwr 



6 



e + 



3(3n + 14)e 



/2 



n + 8 

n + 2 
4(n + 8) 



e + 



(n + 8) 2 

6e' 2 /n + 3 3 
h In - 

n + 8Vn + 8 2 



(18) 



Clearly the difference between the weak-long-range and short-range regimes 
comes from the difference between e and e' as given by Eq.(15). From the 
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work of || we know that the LRI fixed point becomes instable at a = a s . 
The signal of instability appears however at three-loops. Our work shows 
that already at two-loops a new fixed point develops, driving the pure LRI 
model to the intermediate weak-long-range regime. The a dependence of the 
critical exponent 9' in this regime is linear and is shown in the curves 'd' and 
'e' of Figure 3. 

We summarise now our results. We studied the short-time critical be- 
haviour of the Ginzburg-Landau models with LRI in the e-expansion up to 
two-loop order. We observed an initial critical increase for dimensions smaller 
than d c and for the interaction range d/2 < a < d. We obtained the univer- 
sal critical exponents 9 and 9' of the initial slip as functions of d, n, and the 
interaction range parameter a. The limit in which pure LRI is approaching 
the SRI has been also discussed in some detail. 

Acknowledgements: The authors are grateful to H. Luo and B. Zheng 
for fruitful discussions, and thank C. Untch for help in using the computers. 

Appendix The calculation of Zq 

In order to determine the renormalisation constant Zq, we calculate the two- 
point function < s(—q, t)s(q, t') >, with one leg attached to the initial surface 
f = 

< s(-q,t)s(q,0) >= r°dH < s(-q,t)s(q,t') >^ T^(q,t') , 
Jo 

by using the graphs of Figure 4. 

In these diagrams C^ D \t,t') and G p (t,t') are represented by solid lines 
with and without arrows respectively. The small circle means that one time 
argument is set equal to zero. The factor < s(—q,t)s(q,t') >^ denotes the 
contribution to the two-point function coming only from the equilibrium part 
Cp S \t,t'), whereas the residual factor Ti \q,t') is the sum of the amplitudes 
with at least one initial part C^'(t,t'). 

We write the singular part of at the critical point r = in the form 

r§(g = o,£) = h + h + h + h + h (19) 
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where Ij with j = 1, 2, 3, 4, 5 is the contribution of Fig.4(j). These contribu- 
tions are given by 



h = S(t); 



(20) 



6 J (2n) d 
n + 2\ 2 n 



d d p 



j^G p ,(t,t')C { f\t,t') ; 



/ 4 = (Xg) 



x 



2 n + 2 /■* r d d p d d p' 
~6~ Jo t J (27r) d (2n) d 



G p+pl (t,t') 1 2 C^(t,t')C p ^(t,t') + C^(t,t')C^(t,t') 



n + 2 
6 



d d p d d p' 
(27t) d (2tt) 



iC$\t,t)C®(lf,lf) 



(21) 
(22) 




(1) 



(2) 





(3) 






(4) (5) 
Fig. 4 Diagrams contributing to r^(g, t) up to two loops. 



By using the formulae 



dxx^e-^ = (J,-T{u) ; 



t dxx v -\t - xY' 1 
Jo 



r(z/ + /i) 
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valid for Rez/ > and Re^u > 0, it is not difficult to obtain 
n + 2. 



h 
h 



6(7 



-XgK d V{l-e/a){2\t) 



-l+e/a . 



n + 2\ 



X 



) KgK d ) 

T 2 (l + e/a) 1 
r(l + 2e/<r) 2 



; r 2 (f-^) 

2ae 
(2At)- 1+2e / CT ; 



n + 2\ 2 ^, 2 r 2 (f- e /a) 



6 



2ae 



(2Xt) 



-l+2e/a 



(23) 

(24) 
(25) 



By integrating in (21) over t' and over the length of p', one finds 

n + 2 w _ x9 „„ „ , N r r (i) , T {2) 



h=-—^X(gK d ) 2 T(l-2e/a) -I™ + /f + 0(e) (2At)- 1+2 ^ 

0(7 L J 



where 
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(27r) d x CT (e + x) CT (l + a; <T ) ' 



with e the unit vector along the ci-axis. The first integral is easily done with 
the result 



J- A 



K 



r2 r 



^ p(f)r (<l -») 



In order to carry out I\ , one can use the following expansion of l/(e + x) CT : 

OO 

(x 2 + 2x • e + I)-"/ 2 = [max(x, 1)]~ CT £ [min(x, l/x)] n (-l)"</ 2 (x • e) 



n=0 



where x stands for the unit vector of x and c^/ 2 (x ■ e) are Gegenbauer poly- 
nomials. This leads to 



/dx(x + e)-" = j x t 



x < 1 

X > 1 
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(2) 

which can be then used to find the leading contribution to I\ 



/f = - In 2 + 0(e). 
a 



Finally, we get 



-X{gK d fY 1 



6a \ a 



1+ 2 ln2 _^ + 0(f 
e cr 2 



(2At) 



2e/<r-l 



(26) 

The substitution of Eqs. (20), (23-26) in (17) leads to an explicit expression 
for Tiq (q = 0, t) up to terms of order e. 

We renormalise now according to (4) the (bare) quantities entering this 
expression. For the fields s, s and the coupling constant g a one-loop renor- 
malisation will be sufficient. All the necessary information is available in 
Eqs. (5-7). The residual singularity is then removed by requiring [l|] 

Z 1/2 r dte-^Yfhq = 0, t) b = finite for e -> . 
Jo 

Here the subscript b denotes the expression of obtained above in which 
only bare quantities appear. From this condition we compute Z as given by 
Eq.(9). 
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Figure 1: The exponent 9' for n = 1 is plotted versus d. Curves 'a', 'b', 'c', 
and 'd' correspond to a = 1/2, 1, 3/2 and 2 respectively. 
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Figure 2: The exponent 9' for d — 2 is plotted versus l/n. The curves 'a', 
'b', and 'c' correspond to a = 1.05, 3/2, and 2 respectively. 
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Figure 3: The exponent 9' for n = 1 is plotted versus a. The curve 'a' 
corresponds to d — 1. The curves 'b' and 'c' are drawn for a < a s at d = 2 
and d = 3, respectively; the curves 'd' and 'e' are based on Eq.(18) and 
describe the behaviour in the region a > a s . 
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